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The problem of the mutual attraction and joining of two black holes is of importance as both a 
source of gravitational waves and as a testbed of numerical relativity. If the holes start out close 
enough that they are initially surrounded by a common horizon, the problem can be viewed as a 
perturbation of a single black hole. We take initial data due to Misner for close black holes, apply 
perturbation theory and evolve the data with the Zerilli equation. The computed gravitational 
radiation agrees with and extends the results of full numerical computations. 



The collision of two black holes is, in principle, one of the most efficient mechanisms for generation of gravitational 
Ph ' waves. In view of the fact that the LIGO and VIRGO 0] detectors may be detecting events in the coming years, 
. the theoretical determination of possible waveforms has become of great importance. The fact that data may be well 
below the noise level of the detectors may require pattern-matching techniques which require accurate knowledge 
of radiation waveforms. 

The problem of the gravitational radiation generated by colliding black holes is not only of great importance to 
gravitational wave astrophysics, it has also been one of the earliest applications of numerical general relativity. Smarr 
and Eppley more than 15 years ago, computed the radiation waveforms for the axisymmetric problem of two 

holes, starting from rest and falling into each other in a head-on collision. The importance of this problem has 
motivated a recent reconsideration, both numerical and analytical, by Anninos et al. The numerical work is 
difficult, especially when the holes are initially close together. In that case the radiation is dominated by horizon 
processes. In addition, for initially close holes the radiation generated is relatively small and the numerical errors in 
its computation can be particularly troublesome. The purpose of this paper is to provide a method of computing the 
radiated power generated when the holes start off close together. Our method is based on perturbation theory and is 



considerably more economical than a full numerical simulation. It can also be viewed as a benchmark against which 
^Jy numerical codes can be checked. 

For a full numerical computation of the problem, data used are those for the two throats of a momentarily static 
"wormhole," for which an analytic solution was given by Misner [||. The initial data has a parameter fiQ which can 
be adjusted so that the initial conditions correspond to different values of L/M, where L is the initial separation of 
!_h ' the throats, and M is the mass of the spacetime. 

For values of /j,q corresponding to large and moderate starting values of L/M the motion of the "particles," not their 
black hole nature, is crucial to the generation of gravitational radiation. In this case the amount of radiation emitted 
can be understood with a quasi-newtonian approximation that starts with the known radiation for a point mass falling 
into a hole ||. For values above around /iq ~ 2, this quasi-newtonian approximation is in remarkably good agreement 
with the results of numerical relativity. For smaller values of ^iq, however, the approximation seriously overestimates 
the radiation. When the throats start off at small separation their "internal" structure cannot be ignored. 

It is our purpose here to provide an understanding of the opposite limit, the limit of small /iq (and hence small 
L/M). Below we give a simple and attractive analytic result valid in the limit of small separation for the energy 
radiated. A rather complete understanding of the phenomenon is then afforded by this result along with the quasi- 
newtonian estimate, and the bridge between the two given by the numerical relativity computations. In addition, our 
method gives predicted waveforms, and other features of the wave, that can be useful for verifying the accuracy of 
the numerical relativity computations. 

It should be understood that our small-/io result does not have the same robust connection with a simple physical 
picture as the quasi-newtonian approximation has with the mutual in-fall of two holes from large distance. In the 
large distance limits the details of the choice of initial data to represent the individual holes is unimportant; it is their 
particle-like motion toward each other that dominates the radiation. For the initially close limit the Misner data is 
one possible choice for momentarily static initial conditions, but it is singled out by mathematical convenience, not 
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by any claim that it represents the most natural, instantaneously stationary, initial distortion of the participating 
throats. 

In terms of bispherical coordinates fi, 77, 4>, the Misner initial data H take the form 
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where a is a constant with the dimension of length, and 



<PMis 



E 



n— — 00 



-\/cosh(^ + 2n^ ) cos?/ 



(i) 



(2) 



We now change from bispherical to spherical coordinates R, 9, <f> and introduce a Schwarzschild radial coordinate r 
through R = (r 1 / 2 + yjr - 2M) to arrive at a line element of the form 
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with dtt 2 = d9 2 + sin 2 9 d(j) 2 , and with 
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The square root in the summation has the form of the generating function for the Legendre functions P^(cos#), so 
that F can be expressed as 
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Here, and below, we use the notation: 
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Since the Misner geometry satisfies the initial value equations of general relativity, and is momentarily stationary, 
there exists a coordinate choice T, r, 9, <fi such that the initial data generates a 4-geometry at T — of the form 
ds 2 = —dT 2 + (is^jg^j., and for which dg^/dT — at T = 0. One can make a transformation such that the 
4-geometry takes the form 
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For a constant t slice of the Schwarzschild geometry F is unity. The difference between unity and the Misner form of 
F is the extent to which the geometry in (^|) initially deviates from a spherically symmetric black hole. 

We have therefore cast Misner's initial data in Schwarzschild coordinates. We now need to explore the limit in 
which the two black holes are "close." When /io < 1.36 an apparent horizon at R ss 2M surrounds both throats. As 
/io decreases further, the value of L/M decreases, and hence the ratio of L to the horizon radius decreases. As pointed 
out by Smarr the horizon, and the geometry outside it, should then be nearly spherical, and it is only the geometry 
outside the horizon that influences the radiation sent outward to infinity. Linearized perturbation theory should 
therefore give a good description of the generation of radiation (though not of the highly non-spherical geometry 
inside the horizon). 

To give this picture a mathematical realization we note that the coefficients in (^) then have the small- (io limit 
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In the spacetime corresponding to (^|) we can therefore consider e = l/ \ ln/io| to be an expansion parameter. If in 
we keep only the leading term in e for each £ we get 
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We now argue that each multipole term can be treated individually by linearized perturbation theory. This is clearly 
true for the £ — 2 case. To lowest order in e, that is to order e 3 , the field equations contain only terms linear in the 
perturbation (i.e., linear in 1 — F), and those terms are pure £ — 2. Next consider, for example, expansion to order 
e 7 . This would contain £ = 6 terms linearly, but would also contain nonlinear terms, e.g., from the square of the £ = 2 
term (with higher order corrections to (^) included). But the nonlinear contributions can have no £ — 6 part. (To get 
an £ = 6 term requires a cube of an £ = 2 term, or a product of £ = 2 and £ = 4, both of which are higher order in 
e.) This shows that the I = 6 part of the vacuum Einstein equations is linear in 1 — F, and hence can be treated by 
linear perturbation theory. This argument easily generalizes to arbitrary £. 

Though each of the £-poles in ([To]) satisfies the linearized source-less Einstein equations, they are not all the same 
order in e. We take up here only the dominant term at small initial separation, the pure quadrupole £ — 2 term. To 
treat this as a perturbation problem, we use the notation and formalism of Cunningham et al. [Q except that we will 
omit the tildes over variables; equations from that paper will be cited as "CPM." It is important to note that the 
formulas in that paper, based on the work of Moncrief |^|, are gauge invariant, so we need pay no attention to the 
coordinate gauge of (JTc|) . 

From CPM (11-25), (11-26), we find that the only non- vanishing metric perturbation functions are H2 = K = 
g{r)p{jJLo) where (with the notation of (^)) 
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From CPM (11-27), (11-28) we next find 
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whereA = 1 



4^ . The function ip then satisfies the £ = 2 Zerilli equation || 
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where r* = r + 2Mln(^j — l). [It is worth noting that the function Q 2 of CPM (11-28), calculated from our 
perturbations (|l2]), explicitly solves the initial value constraint Q2 = as given in CPM (11-29)]. 

The form of ip given by ([l2])-(|l5|) is now taken as initial data along with the initial condition dip/dt = at t = 0. 
The problem is greatly simplified by the fact that the only /j,q dependence is contained in the multiplicative factor 
p(fio)- Since the initial data are proportional to p(fio), it follows that the evolved waveform -0 is proportional to p((J-o), 
and the radiated power and energy are proportional to [p(p,o)] 2 . It is only necessary, therefore, to do one computation 
of the evolved waveform and radiated power. The p,Q dependence is known at the outset. 

The initial form of ip (with p(no) set to unity) is shown in figure 1. We evolved these data numerically with the 
Zerilli equation Jl4|). The resulting wave form, at r* = 200 (we use units in which 2M — 1), as a function of t, 
is shown in figure 2. It clearly exhibits quasinormal ringing and power-law tails corresponding to a quadrupolar 
perturbation. The values of the quasinormal frequencies ]To| and power-law exponents jll| are in excellent agreement 
with theoretical values. From the evolved data we compute the radiated power, which is given by CPM (111-28), 
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If p(no) is set to unity, this procedure gives a computed energy of 3.07 x 10 6 . The result energy, as a function of fXo, 
is therefore 



This result is displayed in figure 2, where it is compared with the numerical results reported by Anninos et al. ||. It is 
intriguing that the remarkable agreement extends considerably beyond the small-/io region in which our approximation 
is expected to be applicable. 

The general method is suited to a fairly wide variety of initial data. This paper only concentrated on the aesthetically 
elegant Misner data as an example. The method can also be applied to initial data which are known only in numerical 
form and which represent perturbations of a black hole. To do this would require considerably more numerical 
analysis than for the Misner initial data but, when applicable, would be much quicker and less expensive than the 
fully numerical methods for integrating the nonlinear field equations. 
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Faculty Development Program. 
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FIG. 1. The function tp of the Cunningham-Price-Moncrief perturbation scheme for Misner's initial data. The values shown 
are for p(fio) = 1, and for units in which 2M — 1. For the r* coordinate the horizon is at r* = — oo and r* ~ r for large positive 
values. 

FIG. 2. Time evolution of the Misner initial data (with p(fio) ~ 1, 2M = 1), from the point of view of an observer fixed at 
r* = 200. We see the appearance of quasinormal ringing with the predicted period of 8.4. In the inset we display in a log-log 
plot the late time behavior of the field, which clearly exhibits a power-law tail form with exponent —6 as predicted by theory. 

FIG. 3. The solid curve is the prediction for the radiated energy, as a function of fio, based on linearized perturbation theory. 
The black dots correspond to the values of numerical relativity results reported by Anninos et al. [5]. 



Energy/2M = 3.07 x 10" 6 p(/z ) 2 . 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/gr-qc/9402039vl 
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FIG. 1. The function ip of the Cunningham-Price-Moncrief perturbation scheme for Misner's initial data. The function is 
perfectly smooth, although it varies rapidly around r* = 0. In this coordinate the horizon is at r* = — oo and r* ~ r for large 
positive values. 
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FIG. 2. Time evolution of Misner's initial data, from the point of view of a static observer at r* = 200. We see the appearance 
of quasinormal ringing with the predicted period of 8.4. In the inset we display in a log-log plot the late time behavior of the 
field, which clearly exhibits a power law tail form with exponent —6 as predicted by theory. 
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FIG. 3. The solid curve is the prediction of our calculations for the radiated power. The power radiated is composed of 
two factors, one of them a simple analytic function of fi and the other a coefficient determined by numerically evolving the 
perturbation. The coefficient is the same for all values of fi. The black dots correspond to the values of Anninos et al. 



